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ABS SER ACH 


The problem of predicting the motion of a spaca 
vehicle in the vicinity of € Trojan polut in cislunar 
apace is considered, Specificas1y, it is desíred to 
know if the vebiclo will xemein iu this vieinicy, OF 
diverge from the Trojan point. If the equations of motion 
of the vehicle ave linearized at tha Trojan point, the soiu- 
tion of these equaticne inti sax that the vehicle will 
agsune a periodico motion about the ie point. Fundamen- 
taliy, the problem reduces to determining tche validity, or 
cenge of validity of these solutions. 





Various methods are discussed elong with thair agan- 
ciated problems. Further investigetion “Le proposed along 
the lines of ideatifying ad roducing computational exvoxs. 
Also it is proposed to invest igate che nigniffc sace o£ 
linesrizing the biet ions of motion at points displaced 
from the Trojan pointe 
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OBJECT 


It is desired te study tha orbital stability of a 
Space vehicle et a Trojan point in cislumar space, 1.e., 
waere the earth, moon and vehicie form en equilateral 
triangie in the orbital pisme of tha eerth and meon., ££ 
@ vehicle originally placed at tke Trojan soint is indeed 
in a stable orbit relativo to the Trojer point, it should 
then be determined whet order of inaccuracies of injec- 


tion can be tolereted to stil: establish a stable orbit. 


This aasumes that there is an crerational requive- 
ment for a space vehicle of this sort which maintains ite 
relative position with respect to the earth and meson. The 
reguírenent might be for e space navigetion or ecumunice- 
tions setollito, or it might be for a replenishment station 
for e manned deep spece probe. However, the specific 
advantages of e Trolan peint orbit were not investigated 


further. 
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cipliunar vohicle with lees than escape velocity will remain 
in eislımsz space, barging collisions with auntkar body, end 
wili kave an upper bound to ite displacement fron suy point 
in the earth-moon coordinate system.  & comvenient rostric- 
tion on the rigorous definition of stability le not immediately 
apparent. However, the dietinction is perhaps overly meticu- 
loug. We sre concerned with s vehicle whieh remaina in the 
vicinity of a Trojan point and, due to initiel displacement 
from the Trojem point, initial smon-zero relative velocity 
with respect to the Trojan point, or pertubations exterral 

to the system, appeara to asclliate about tha Trojen point. 
If it centinuce te do this, it i8 Stsblo; i£ it diverges from 


tho Trojan point, it ia unstable. 


Historically, the Trojan polacs save been considered sta- 
ble for "small disturbances” based on a iimearigation of the 
equations of motiom of the thira body at the Trojan point. 
However, the etobility of a body im che actual nom-limear systes 
has not been comelusively determined. In feet, there are 
arguments indicating thet, oven granted exact knowledge of the 
physical constente ¿nd unlimited precision and spesd of digital 
computetion, the stebillty can still not be established with 


mathematical rigor. 


The existence of the Trojan asteroida is not as good an 
argumant for stability as it seems et firet thought. Aeccordiag 
to the solut£on of their linearized equations of motion, the 


periods of their two modes of vibzation ere 11.89 yesre and 





407 years. Te corresponding portlods in eislumar space ere 
28.6 days and 91.7 days. Thus the 54 years evailable for 
observation of tho Trojan asteroids correspenda, in one 
sense, to only about ii days observation of a ciglumar 


Trojen point satellite 


We would like te knew specifically the ranges of initial 
conditions end magnitudes of perturbatione for whieh a Trojan 
point vobicle would be stable. However, te study described 
here does net reveal such & eriterion. Pather, it seems we 
muet be content with specifying e time Limit during which the 
vehicle will aot exceed & spccifiod displacement from tho Tro- 
jan point. The imtent of this invostigetion was to provide 
&nswers Co questions of tha sort, What are the guidance 
accurecy reguiremsats for establishing a cislumar Trojen point 
space station to remain within 29,00) miles of the Trojae point 
for et least bo years?" Again, it must be admitted that the 
&nswexr i5 not coceinse iu the present investigation. However, 


the diractions toward whieh such answers ile may ba indicated. 


$ 





CHAPTER 2 


RESTRICPIONS ON THE PROBLEM 


The problem ie considered to be that of the restricted 
three-body problem, where the mass of the third bedy, the 
vebieio, 15 negligible compared with the other teo bodies. 
The vehicie will thon have no e££act on the motion of tha 
other two bodies. Thie i5 obviousiy valid for the problen 


of earth, moon and a vehicle leunched from the earth. 


Another restriction is the elimination of all portur- 
bative forces extermal to the esrth-meon ayate. This is 
equivalent to considering tha berycenter of the carth-moon 
system as an inertial point, end neglecting the curvature 
and velocity variatious of ite crbital motion. The partur- 
bative effect of the sun is perhaps sesential in che eolu- 
tion of the physical problem. However, it can later be 
added to the differential equetions es a forcing function 
when the Simpler system le more £ully undorgtood. Other 
perturpative forces such es that due to Jupiter, tidal 
forces, earth oblateness and radiation pressure ara probably 


indeed negligible. 


¿OXe 


Ley 


Aiso tha orbir of the mon about che carth, or, 


un 





precisely, etour the berycenter of tas carth-290n System, 
i5 congidered to be circular. Tris effect can also be 
added leter 1£ desired. 
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Erzors duo to those siz 
álecuesed further in the text. Their effect is also cen- 


giderad by Buckeim (1). 





CHAPTER 3 


ee. 


ERBE LIERATLON POXSTS 


2.1 The Geometey of a Epace Y Jehiele , &ttrocted bv Two Bodies 


In Fig. 3.1, the vebicle, Mgs ie attracted by the two 
celestial bodies, My and m 2» mich revolve about their inar- 
tial bexvcentor, B. Tn eigluner space, Ey and g. axe tho 
earth end moon, and, Zur comvenlenes, they will bo sc namen. 
The mass of tho vehisle te considexed small enough te nave 
mogiigible effect on the motion of the sarth and moon. 


+ 


To simplify tho cqu&tiong, ell distances ere normalizges 
wich regpect to the esrth-soonu distances. Distances Erom 
earth sud moon to their baerxycenter are thon © and A regpece 


tively where 


The grevitational forcas on the vehicle sre fenreseutes 
2 


7 E + 7 tpa D» e g” a ni Ese = . ES ton 
wy £, and £o which are vector specific foress, ox forces per 
c» 


vnit messa of tae vohicie. 
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lt ís convenient to choose a rotating coordinate 
eyaton la walch the earth and moon appear stationary. 
The equations of motion in this ccordinete system must 
then include centripetal and Coriolis eccalerations 
resulting from the coordinate trensiczmation. It is also 
convenient to consider the centripetal acceieration of the 
vehicle due to coordinate transformation as a centrifugal 
specific force, È o: along the radius vector from the bary- 


center. | . 
£, = x" 


where ía the engular velocity of the esordinate system 
with respect to inercial space. Then from Kepler's third 
Law 


Q~ s 


ye 


(3.3) 


fo 6 (adm) x 
3.2. The Existence of Libration Pointa 
The following development, due to Klemperer and 
Benedikt, (2) shows the existence of five points in the 
plene of the moon's orbit eb which a vehicle with zero 


relative velocity in the rotating coordinate system 


experiences sero resultent force. These pointe were 


DB 





first discovered by Laerenmgs in 1772 as particular 
solutions of the three-body preble::. Their existence 
could be proved with lesg effort by simple substitu- 
tion of their coordinates into the equations of motion 
but the Following development also shows that these 


are the only five pointe at which this condition cceurs. 
Consider the vehicle to be at rest in the rotating 


system, With zero relative velocity, the Coriolis 


acceleration is zero end, 1£ the vehicle is in egulli- 
brium 
f, +i, +4, 20 
Consider £ivat the circumterentiel equilibrium of 
the vehicle, or the sumeation of foree components par- 
pendieculer to tho radius vector from the barycemter. 


~£, sin y + f, sinye g 








AG (mpm) EG (my) — 
- pá aia t —— sin Yoo (3.4) 
r Zo 
A pinya =p sin y 
aot al 
> Y "5 
One iocus o£ ciremferential equilibrium is the straight 


lime through y end Ein RETO 
sin y e oin Y = G 


Also, by the Law of gines 


x9 





Bingo E Sin i ( p". 2) a8 ixi e 





Es CE 
sin Y = Bin $ 
a 2 


Substituting for ain Zand sin fA in Eq. 3.6 


ALE sin? ) 3 € Gain $ ) 
3 3 
r $ 


fad the other iccue of circumferential cquili» 


Ed e 


Trium is 
the line 
Ea = Eo 
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Hext consider the radial equilibriun of the vehicle, 
or the timeation o£ force components peraliel te the radius 


vector frem tha baerycenter. 


-£4 cos - £, cos M+ f° 
tg uu tis 
B x61 2 cosg Eus “y an hm COB 2 + TG (m4. ) e 





zo v4 
u 2 
t cos 7 ~ € cos?’ Prec 
d A f 
j “2 
y (3.5) 
Finding cos Zand cos by the las of cosines 
el o gl d r$ «> ZEE ecs Æ 


cos Za qian By” Ex 


qe T ee " ZEE, COS ff 





ı 2 M =, we um, 


Muitiplying by -26z7 => to clear tho fractions 


2:3 3 


. US 
ET HAEE, „Re? + Er rj * ex, 4 x 


x -eR r3 -2r Eiko 9 J 
(3.6) 
To Gliminare r from the equotion, r i5 found in terme 


o£ ty and Es by the law of cosines. 


2 A E 
e^ 4 xT e 2er, cos 2 (3-2) 


and alge 

x - l4 r? - 2r, cos (f (3.8) 
Substituting for, ces “yf from Eq. 3.7 into Eq. 3.8 
p? = e? 4 zt + E (iter) set + c (1-8) + cas 


v? » e(1.6) & x2 (1-0) + er? 9 6r AT + ext 


Substituting Lor p? in Eq. 3.6 end combining terms the lecue 


for radial ewuilibriva is obtained 


E AE 2 3 
Hidra MPO 
+27 > v2 a p 


The equilibrium loci are shown im Fig. 3.2. The five 





intersectiona of the two sets of curves are the only pointe 
4 


at which e body would *xceuoin at equilibrio. These five 


system, but the development; ylelde no dyaemic information. 
To specify tne behavior of a vehicle im the vicinity of a 
ALibretion point with deviations im position end non-zero 

relativo velocity, the equations of motion of the vehicle 


must be investigated. 
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Fig. 3.2. Location of the Libration Points 





CHAPTER 4 


THE EQUATIONS OF MOTION 


4.1 Formulation cf the Equations of Motion 


In en inertial coordinate system, the motion of the 
vehicle can be described by eguating fte acceleration te 
4¢3 gravitational esecifie force. Sinsa 26 Lea desired te 
describo the metion in the creriiinké wet whieh rotetes with 
the moon, the additional apparent accelerations which appear 
in tho coordinate transtoxmmation must be included. Tha total 


acceleration congists cf the gum of the following teres: 


L. The aceeleration o£ the coordínate system in 
inertial space. This is considered vero in the restricted 
three-body problem. In the phys£col problem im cisiunar space, 
it corresponds to the seceleration of the sarti-mccn system in 
ite orbit around the sun, including the perturbation effects 
of other planets. Ets only significant component is tha een- 
tripetal acceleration of the earth valeh seems to be the most 
appreciable effect noglostod in confining tha investigation to 
the restricted theeea-bedy problem. If later it appears desireble 
to congider this term, te eaz bo Included as a perledic forcing 
function im the auuntions of motion, representing tha differenes 
between the specific aravítationel forse of the sum ou the 


TA 


D^ ! 





vehicle and that force on the baryeenter of che earth-meon 
Byte. 


2. The tangential ecceleration. This is also 
considered to be zero since the moon is assumed to be in e 
uniform circular orbit about the earth. If tha effect of 
the eccentricity of the moon's orbit were to be considered, 


this term would kave to be included. 


3. The centripetal acceleration. This isa treated 
an a centrifugal specifie foree and grouped with the gravite- 
tional forces. 


6. The velative acceleration. Im the (x, y) rota- 


ting coordinate system, its components ave E end y which ara 


included in tha equations of moticn. 


a 3. Tre Coriolis acceleration. Its magnitude is 
at £g 4 y^, and ite divection is defined by tho vector 
creas product of the system angular velecity and the relative 
velocity of the vehicie. Its z Bud y components are -299 end 


TUE respectively. 
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ás ahora in Fig. 4.1, the origin of the csorsdinate 
system is chosen at en equilateral triemgle point, oF 
Trojan point, 5ince wo are primarily concerned with the 


dyneanies of & vebicie in the vicinity of this point. 


The x and y componente cf the vecter equation of 


motion are thon 


A 
Substituting from egustiens 3.1, 3.2, and 3.3 and reaolving 


> 


x - 20$ = 26 Qm, ) QUOS y + i» Se co B + Gn tm.) Y cos d 


a 


z. ^ 
” Kae) ann Com x i ile 
y + 20% a ACM) gin ys SC ee sin B + Giyim) z nim 
2” x 
1 a 


To Zurther sispli£y the equations, time will be noxmealives 
with respect to the period of rotation of the earth-mocn syscem 
so that: 


£3. «los, m 
el e GG.) (2 were t' is now real cime 


ierat 


a? 


" dx "uh. w dE 
> 2f a cf. o T 
ten gE ~ Verte) g 
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E q ge mE, = Slate) La, 
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with the dot demoting diiferentiabien with respect te 


dimensionless time. 


In dimensionless time, the equations of motion 
become 


x > 2 y ES L cos Y4——, cos f+ Y Cog P 4.2) 
| z 


r 
2 


29 fi 
y +2 TETE A sin yy - - Sy aint z oin $ 


1 E (6.2) 
Referring to Fig. 4.1, the quantities Eq Tys and the 
trigonometric functions Of? » P 3 end D may be expragsed 


in terns of x end y as follows: 
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24 2 E 
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y 3 ER um 2 
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' T pe 
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Substituting in equations ¢.1 and 4.2 the equations 


of motion become 


-3 a Y =, 
k-29 = xeh-6- (sets) (y^ deca. (:t-235) (y^ Jy a +1) / 


(4.3; 

T, E A -J V4. We -Y, 

y 2x € m - (yt Pt (y^ yks eil) / s (yi > ye) z 
(4.4) 


4.2 Integration of the Eyuations of Motion 


* 


The only known integral of these equations was discovered 
by Jacobi. It has since been proved that, in this coordinate 


system, no other algebraic integrals exist. 


T ^ 2 
Let v = -Cx SG oye el 8 21 GI ye ee) 


= 
T (y^ yes? etl) l 
(4.5) 


yg” - / 
Thon È e -erot to nd) /2 


Ox zo 
-€ (x-L ) (ya V pat +1) 2 
(4.6) 
| ne X 
dl a “A or 3) (y^ a yix firl) 2 
dy -3/ 
-E (yF 26 + Sybex? ash) 72 
(4.7) 
and the equations of motion are 
x-2ye Qu (4.8) 


Ox 
LE 


n 





^ m è C 
m m e 
3 vo P 32 TZ ugs (e 2 $) 


Multiplying equations 4.8 anc 4.9 by 2x and Zy respectively 
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2xx-2yxmo2x 95 (6.10) 
OR 
2'5$ 2yx 2 25 29 (4.11) 
: dy 
Adding equations 6.10 sad 6.11 
zx + 2 y y. 2 £8 + 25 n 
Ox . Oy 
2x SÉ 4 25 49€ 2 g 22 dx, ¿22 ay 
de . “de ex dc. sy de 
Integrating 
RPE 222 | ¢ Beas Say) = zu+ 6 
dx . OF 
sa 12 
Ge «2u d m Y y 
The Jacob! integral is 
Hu 
G e 26x - (y y^ m ibi )-2a(y yia 1) 


afk Cy s E ey el) fs $° 2 J^ c COREL 
| (4.12) 
6.3  Vestorial Espresentation of the Equations of Mosion 
7: Should be noted here that the subsequent: 
equations apply to oniy one of the two Trojan points. For 
the other Trojan peint, e syametry sbout the g-asia applies, 
bue the equations and directions may not be applied ilteraliy 
without sonsidoraticn of Chis symmetry 
EE the motion of the vehicle 23 represeated in che 


comples plene where 


A 
K 





E 
g 
n 
+ 
pu 


equations 4.3 and 4.4 con vo written s6 


x-2y "£f (4.33) 
ly + 2ix = i£, (^, 34] 


Adding eqvations 4.13 ad 4.14 


2 + Zig = f+ iE, (5.13) 


At first thought, i6 would geen thet, for a stabile 
Trojan point oroit, the resultant of Ez ad E, ahouid be 
directed toward the origin. This is mot the cese; i6 will 
be seen later that, in tho vicinity cf ehe Trojan point, £ 
aad £ i. kave the gene siga as x end y respectively. Rather, 
it is the GCorielia acceleretion whien provides sa average 
velative acecleratica, z, toward the origin if the orbit is 
to be static, The vectorial reoresontetion of Eg. 4.15 is 
ahown iu Fig. 4.2. Thie indicates thak the orbit muet have 
@ gross motion cloekwiea about the origin fur stability. 


Au A n 23, 
i s wer 


SE ene ee ee 
<r 
KM 
^y 
¿o 
DA 
£e 


ro - è 


PATEAR IDE CESAR 
$ TERRENO TORES TETAS A Tinte Nip TA RR AMT Oe X 


ü 
Fig. 4.2. Vectorial Representacion ol tas preciosos of Koki 
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4.4 Linesrizastion of the Equations of Motion 


As a first order approximation of the bonavior o£ & 
vehicle ia che vicinity of a líbration point, the equations 
cf motion are linearized and the solution of the linearized 
equations ig examined. The procedure used by Moulton (3) 
is essentialiy to linearize the right hend sides of equetions 
4.1 and 4.2 by Taylor series and solve tbo resulting linearized 
equations. Whittaker (4) uses a different appreach, working 
with the equations of motien in the form of four first order 


Gifferentilal ecuations, with the seme results. 


The first approach wifl be extended kere to obtain a 
more complete linesr model of the system, including an integral 
of the equations equivalent to the Jaccbíen. This will later 
be used for e more complete comparison of the actual and 


linearized systems. 


Linearizing the differential eguetíons i58 equivalent tc 
retalaing second order tezms in the Taylor series vrpre- 
sentation of the Jacobian function, Eq. 4.5. Letting primed 


symbols represent ths quantities ef the Linearized systoa 


Eq. &.5 can be written as 
2 


y 2^ Qu ci £ à ce e 
u' (x,y) 2 v (0,0) 4 — (0,0)z- <= (0,0) ye; = 5 (0,0)x 
eut OF à Oz e 
(6,16) 
2 2 
„U {0 ,O) xyes IE (o, 0)y* 
Judy. aye 


Substituting z = y = 0 im equations 4.5, 4.6 and 4.7 
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u (0,0) = ) +74 € "T 
2 


-$ (0,0)= -€ + 3 - kre bc eg 


88 (0,0) « '3 .3 v3 _ DET 
2 Ra 2 


Differentiating equations 4.6 and 4,7 


2, A -3 
= E z) yyte ttl) Y -AG^* Byra o eti) /2 
% 


ar 
£36 (x-5) ^ (y^ Gyan” -X m) "v eG (y + yr -xbi) EZ 
(4.17) 
2 = 9) y 54 NC 
=a = E. ye By) Saly e Byta ttt) y 
y P mA E mé -Fo 
23 (ys 3 y apta xvi) /£e Cy ep 1) "A 


(4.18) 


2u e E 
E æ JRL) (yt & (y^ d ut ah) “hs 3€ (20-25) of (yA ¿By? Župy 
xOy 

(4.19) 


Substituting x = y = 0 ín equations 4,17, 4.18 and 4,19 


& . ; : 
u (0,0) a % 4 2 A - Rr 3 C-G um i 
ax? ly é 
un (Q Q) e1422-242 6 52 

2 = = wil 

Ay: e. DEG és 

wu. wn » NS 
Su. 35. a. 363 (3.52) 23 (1-20) 
dxdy & 4 4 ty 


Substituting thear quantities in equation 4.16,the Jacobian 


function for the l£ínonsrized model le 





u’ e + : x^ 4 aa £1 20) xy + y (4.20) 


tu Ju 


Differentiating, the linearized equations of moticn are 
ootained 


ng . du' 3 343 
z -2y e Mn 2x4 22 (1-28) 
n Ox & & á 


(4.21) 


o o > op! 343 (1-28) x 9 y 
$ 


(4.22) 


fad the interral of the equations, equivalent to the Jacoblen, 


£8 


«9 «5 e V3 2 
= + 


(4.23) 





ENE SESNTTOS OF THE LIREARIZRD ECUATLÍONS OF MOTION 


2.4 iva Stability Criterion 


in differential operater form, equations 4.21 
and 4.22 can be wrítten as 


& n & = 
y 
[2p - 313 (1-20)| 2 (p^ - 2) y 2 0 
és 4 
The characteristic equation of this syeten is 


Ya aya 
(p^ > 2) (p* + 2) + E + EM: (1-2) [2 . xd avail Q 
5 és & - & 


oZ 


The voors of che eharactor£ictio equations ara 





2 censi Herr 
PSA LS e fá 
$ 


Zé the second term te veal, the crerecterístic squation 

yields two pairs of puro dmasinesy Looks, end Ghs solution 
consists ct pura @ismgoido. However, ££ Che second term ie 
imaginavy, tne vooÉs kave positiva zesl parta, and the ame 


piitués of the solution ineresses exponentially with time. 


23 





The stability criterion is than 


0% 


| 
in the earth 


9: a $ (1-6) 70 


2 1 
faar” Cmax * 23 79 


& 2a eN dz « A, 2 0.0 
a * Je = 0.03852 
fy \ 
f f 32:7 n ls 2d e 96 
inex 


emoon ayetemn, E is approximately Sis and 


the stebillty criterion ia mek. 


5.2 The Straight-line Libration Points 


Movlten used 


a similar devslopment to show that the 


| Straight-line Lisretion pointa ave unstable. Another presan- 
katlon of this approach is made by Buchein (5). 


If the equations of motion era lineerized about a straight- 


dine libration point, the enaractesietic equation of the resuit- 


ing system yields 


two Amaginary roots end tuo reat root8, 


one of which £s always positive. A park of ts solusion 


must then ireraasa expenentially wish timo. Tha system 


is thus considored wnsteble; a swell disturbance in displacsmene 


fron the idbration poimt will imereage with time. 
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3.9 The Trojan Polne Solution 


Kf the mags vatie lo guall enough, en $t io in the 





carth-aooz systen, the linearlea!d equations of motien yield 
@ Solution of The form 


(5.1) 
ya be, 82.2 Pot + Bo 208 Pat + Go 221g Pot $ D, COS pot 
| (5.2) 
where 


2 ARSS DO CUAL DA ZA exit 


pi keki-27 (1-6) 


2 MARITIS EF T -EY eR EESTI “Y in 


po" b kYi- 27€ (1-6) 
aná four of the eight constants of interration are independent. 
The reletlonsbips cong the constants can ba determined by sub- 
Stituting the solution back Into tha linear ditfereatiel eque- 
tions. Solving for the constante of Eq. 3.2 in terug of the 
constants o£ Eq. 5.1 


" $3 (eze) A, - 2p,E 


Ag $13 nen ml e a (5.3) 
e 0*2 
e 


9 
$ 
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538, 7 2 13 (ue2*) B 


t 
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be meins Ban u eR tn a e AT ee Se 
E ne wpa =o 
4 à 
E . 
hinge a 
Gu cmo 9 4 
> wry ¡veran AIC ¡DONA NA We ww Um CA o en 


P3 


2D4 E 


L = 3 (1-28) B. 


y 2 S23 CAMERA SR Ee RD OTN "amate NEGO 6 ago Eisen A M REDE Oa S 


ae 


tv 


$ 


pe + 


Á 


The consterts À., B., 


x» 
* 


orcitiong e£ dfopls 
y «9. 


O., D, are Chan 


the inítial c 
2] (9), ( Q) 


^ 
-* 


e 


8.5. Bliuminetios o£ Que Mode nf Osclilerion 

To bencar umlerstand tip nature ok the 
has been ebtaioc?! Zcom the linearinad ugentions, it 
be detizeble to Jsternine hm the solution is divide 


the two modes of vibraticn sf a Engutioa ef init£a 
Tae constant ot equeriono 5.1 amd 5.2 em be 
a funecion of initial condicione, im tenons o£ Los 
fregueneihas D. eud Poo bis ig dese in Appendix A 
mmapies) coustencem For tbe eaxth-rconm system, but 
physical prostat by biek tke syoten ls excited to 
proporticss of tbe two modes of vibretion 12 net rad 
appavcnc. 


% 


$ 


acement aed velocity, x 


(5 35.6) 


er 
un 
1 


atexminsd by 


(0), 


al condir iom. 
Gotexmised ps 


natural 


USAS 
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To obtain further insight along this line, we determine 
what specification of initial cenditicns will excite only 
one modo of vibration, say tho Pz mode. Foz gimplicity we 


arbitrarily specify the imitial displecement as 


y (0) = 0 (5.8) 


We are, of courge, concerneó with auch smaller displacements, 
but the resuits can bs scaled down. Lf wa epecify that the 
initial velocity ia chosen to eliminate the Do mode, then 


equations 5.1 and 5.2 become 


x = fy Sin p.t + 84 cos pyt 

y e Ay sin p,t + By cos pyt 
Substituting inittel conditions 

X (9) » B, =Q 

y (9) = Be l 

x (0) ® Pty 

y (0) * pyA) 


Also from Eq. 5.3 


ar 


£o 
iur 
1 
e 
pa 
4 
wD 
m 
nec 
> 
qu 3 


fg = 


"u$ 
wh RD A tO 
«> 
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B, - 8; 


Du 





2 
P1*7 


Since B = 1 then A, > 1, (p? 4 2) 
294 & 


Thus the initial velocities must be speci iod as 


& (0) e 5 (p^ 4 2) (5.9) 
m a 

y (0). - È 

y . 3 (1-20) i31.) 


And the eubsequent motion ef the particle ig specified as 


Z 
xo Py + sin p,t (5.41) 





3 
en \ 
4 Ya 1-29. gin pyt + cos pot (5.12) 


y ta 
2p4 


Solving for sin pyt and cog p,t in terms of x and y, we 


can then eliminate time by the equation 
Bin ga e 2 5 a 
[8 pyE J (x,y) + cos” pt (x,y) œ: 


After simplificatlon, the equation of the "orbit" about 
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the Trojan point is found to be 


y 


L3 2 

e - 2 

2 + 3 3 41:22) xy + 27 (1-28) + “pa „2 ad 
2 (py + 2 1 


16(pÍ 4 2 ) 
4 , 
(5.13) 


This is the equaticn of an cilipse with its center at the 
origin. It shoul’ be noted that this is not a Keplerian 
elliptical orbit. Anguler momentum ig not conserved in 
this coordinate syatem, and there is no central force 
field analogous (o that oX a Keplerian orbit in inertial 


spece. 


Although th: motion of the vehicle about the Trojan 
point described in the relative coordinate system may not 
be correctly terred an orbit, the term will be retained 


fox convenience. 


If we let v be the magnitude of the velocity of the 


vehicle in the rulative coordinate system, i.e. 


2 s? oy 


Then this velocity can be found at any point on the orbit 


from Eq. 4.19 


ve a + 3 + x^ IE (1-26) xy A» y" 
t 
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C' is a constant vhich can Be caleulated from initial con- 


ditions 


a 
gn 


Cl y* tu) - 3 


From equat£ons 5.9 and 5.10 
1 


2 
v? (0) «X (d & 2 +27 (1-2e)” | 
64 
ana at emy point on the orbit 
2 m. 2 2, 27 9 dae 
y ez (py +9) +4 (1-26) -— Hein 1-96 
E E 4 Tg aa ae m 
o » 
i 4 ye (5.14) 


Considering an arbitrary point on the orbit to be the 


^ 


initial poiut, the initial velocity necesgary to excite only 
the p, mode i8 determined., Ita maguitude is cbtained from 
Eq. 5.14, and ite direction is specified by the tangent to 


the orbit computed from Eq. 5.13. 


LE we now divide ail displacement and velocity varla- 
bles in equations $.7 “bronsi 3.14 by an arbitrary scaling 
factor, we obtain a new orbit with the same center, shape 
and orientation, which is alse a Soiutioo ot tbc linearized 
egustimaa with imitiel conditions specified to eliminate 


a 
the Be mode. For a continuous varietion of ecaling factor, 
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we have au infinite femíly o£ slileses covering Ihe x-y plsze. 


z(0) em ) Q (3,25 

yO) =k (5.16) 

RO) - > (py = (5.17) 

$(0) » s 2K 3, (1-2€) (5.18) 
pe q 2 

y= E TE 3 (1-26) aim p, + 5 cos pyt (8,20) 





y 
y^i d rc ja t petat «Js iE 
aet i 2) a + 2) k 


o 
t} 2 & yj m 2 2 " 
v" eu nf pe o 2) ehe 27 (1226) e EL: eze i M ed > 
& d . 645 & & - 


es & 59 
a (3(1-20)5p+ E e 


> 


Given any yolut in tho x, y piano a8 tke initisl poipt, 





we can caiculate k by substitution ia Eg. 5.21. We can 
then determine the required initial velocity vector to 
excite the Py mode only from equations 5.17 and 5.18 and 
the resulting orbit Le specified by equations 5.29 through 
5.22. 

If we wish te eliminate the p. mode instead o£ the 
po mode, the procedure is che game, and the resulting 
equations are identical with equations 5.7 throush 5.23 
with the substitution o£ Py for p, 


3.2 Determination of the Orbis £eom Specified Initial 





Conditions 

We now have the mathematical equipment for answering 
the question, "Given an arbitrary initial posítica, what 
vector initial velocity ls required to excite only mode 
p, Ox py?” We do mot have as convenient a physical repre- 
sentation fox answsriag the question "Given an arbitrary 
initial position and an arbitrary initial velocity, what 
does the orbit look likev" Tha texptation is to substi- 
tute rumbers into the equations of Appendix åA. However, 
we can do it sewi-graphically by e@ two etage trial end error 
method, aná perhaps gein a better physical iucight of the 


problen. 


Let thea esuetions of the orbit be written as 
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x= x, + Xo (5.23) 


where x” Ar sin pyt + B, cos p,t 
y; * À, ein p,t * B, cos p.t 
x, " C, sin Pot + D4 cos Pot 


To determine the orbit resulting from specified initial 
conditions, pick a Soie [x, (0), y, (€ as a trial. This 
determines an ellipse of the family of Eq. 5.21 and specifies 
the required x, (0) and ¥, (0). ft also specifies the point 
[ x, (0), y, (0)] by equations 5.23 end 5.24. In turn, this 
point also specifies an ellipse of the family of Eq. 5.21 
with Po substitut.d for pi; end gives the required k, (©) and 
740). If the trial is correct, the vector addition of v4 (0) 
and vo) gives the specified initial condition v(5), and 
the orbit can be constructed by the addition of a vector 
generating the po ellipse at a relative angular velocity o£ 
Ea 
P1 

This method i5 illustrated iu Fig. 5-X which is not 
drawn to scale. It i9 included only for possible clarifica- 


tion of the above explaneticn. 
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5.0 Desirable Ranges of Initial Conditions 


From the preceding amaiysio, it appears that, given an 
unavoidable displacement error, wa will ebtain the "smoothesc" 
orbit if the initial velocity has a magnitude in a renge auch 


that the -Jacobian) constant o£ the vehicle is approximately the 


> ? 


5 


same 83 thet of a body at rest ot the Trojan peint. Ic can 
be seen by substitution in Eq. 4.19 that a body vibrating in 
the p, mode oniy has a relative energy higher than that of a 
body at rest at the Trojan point, while a body vivretinag in 
the p, mode only’ has a zela tive energy lower than that of a 


body at rest at the Trojan point’. 


More inportant, we wish the initial velocity to be in a 
direction compatible with & clockwise orbit about the Trojen 


point (using, the notation of Fig. 4.1). 


t=} 


he method of the preceding section shows that, for ini- 
tial velocities in the “eorreet” direction, we tend to get 
contributions frem the two modes which are additive, and the 
outer bound of the resulting orbital motien is not much different 
Erom the initial displacement. Gn the other hand, if che ini- 
tial velocity is in the opposite direction, ue tend to get 


large contributions Gf the two modes which subtract st the ini- 


tial point, but add at the outer bound o£ che orbit. 


This veri£ies the impression we ger from section 4.3. 
Since the a and centrifugal force field is, in general, 


directed away from the Trojan point, the Coríolis acceleration 


^^ 
$ 





must provide the net neceleration toward the Trojan point 
for a stable orbit. Ver this to eceur. the gross motion of 


the vehicle must be clockwise as can be seen in Fig. 4-2. 


Fig. 5-2, whict is also not drawn to scale, is a sketch 
o£ two attempted in ections of a vehicle ints a Trojan point 
orbit. in each caer the approach is from the negative x 
direction, and she vehicle has been provided with approximately 
the correct relative energy. The closest point of approech 
to the Trojan point. Is considered the initial point for 


generating thu orbit. 
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Fig. 5.2. %lece of Diepldcement Error in Establishing a 


Trojan Point Orbic 





Under these conditions y a miss to the left of the 
Txojan point results in an orbit whose maximum displacement 
from the Trojan poiut is much smaller than that of the orbit 
resulting from a miss to the right cf the Trojan point, In 
an actual attempted injection of a vehicle into this type of 
orb:t it appears we should evaluate the expected angular 
error in guidance , and then aim to the side which will insure 


& reeulting orbit of the first type, 
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COMPARISON OF THE ACTUAL EQUATIONS OF MOTION 


WIIN THE LINZARIZED EGUATTONS 


A 


We will hete, E. he residual force, to be the 
veecor difference betwean the actual feree exerted su 8 
vehicle by the combined gr: vitational and centrifugal force, 
aná the fictiítioue force wi ich would be exerted on the 


vehicle by the lincsriced podel o£ the syaten. 
+ 


at - T e guadu- ged u’ = grad u 


E 


where 


u 2 yy! 


and che quentíties u and u’ are defined by equaticna 4.5 
end 6.6. 
WA 


e$ 
Bartual 
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Figure 6.1. Ef£ccr of the Residual Forse on the Acceleration 
of thy, Vehicle 
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2 TA un me v ^ ^ : ; , ¿ve > E 
Tho xohlduci £orzoo dizectly roeorescats the effect of 


E 


equations cf motion wick 


en 
et 
a 
i2 


ail higher order teme o 


were nogLi:cted iu ^hoiz linosrizution. As ghown, in Fig.6.1, 


v 


if the vahicle is agaumed ta foliow the solution of the 
linearised equations ag a firet approximation, then the resi- 
duel force vector shows the tendency of the actual systen 


to cause Che vehicle ca deviste from thie path. 


it would be convenient Lf the residual force field 
were found to have come significant pattern. Fox exemple, 
££ a goluticn of the linearized equations was formed where 
the residual feres was always inward cawamtithe Trojan polar, 
thia would ba an indication of the selcience of a stable 


orbit somewhere within this regica. 


Several points of the residual Jacobian function, u, 
were calcuioted by the methods shown ía áppendiz B. The 
resulto of these calevietiens ave chova in Fig. 5.1. The 
residual force is the gradient of the function a, $i.» 
it acts along the dircetion of grestest rate of ineresso ox 


‚to that rate, of 


w . e 
u with dísteuce, and ita magnitudes ig equa 


x) 


f 


inereszse. The oiet of Fig. Bed is neither aurficleatly 

complete nor accurate enough for a numerical determination 
of the residual force. Im any cese, it ig more easily caleu- 
lated divcetly by component: with fhe methode of Appendix B. 
instead the plot wes meant to give a reugh indication of the 


direccion and meznitude of the regiduai force es a function 


of displesement Freon the Trojan point. 
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Ho Significant pattern of residual force distribution 
was found. If we trace out an orbit about the Trojan point, 
we sind that the residual force flunctuates in magnitude 
and direction with no readily observable net result. We 
also sea that the closer we are to the Trojan point, the 


smaller is the residual force. 


So far, we have found nothing to indicate that our 
linearized model is not a good approximation of the actual 
System. The deviations show no significent patterns, and 
they seem to decrease regularly as we get closer to the 
point about which the equations were linearized. However, 
we are still not gumantésd that an orbit about the Trojan 
point is absolutely stable for "small disturbances", nor do 
we have any information about how small the disturbances 
must be to obtain "partial stability,” i.e., keeping a 
vehicle within specified limits of displacement from the 
Trojan point for a epecífled period of timc. 
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ESAPIVA 7 
AER EnUDLUM OP DUIUESMINIUG DEADILIIV BY DIGITAL 
COBPUITAXION OF Ths OSSYY 


Usually vien we seal an an@lyticeal aslution cf a system 


oZ non-linosr diffe 


ch 


contisal ecuetiona, we feei that wa can 
aiwaye fall back on numerical method is all else fails. We 


would like to have a scluziom im closed form or, in this cage, 


a 


an anelyeciícal upros bend on che orbit, but, 23 general, this 


y $ 


. 148 & convenience which saves o the Lins and expense of resorting 


to a digitel computer. Here, khere 1a some doubt about whether 
we can obtain a saclelcoloz; amewver to the stability cuestion 
t Au 3 
by nomerical mothods. 
Tuc obvious methos is merely to inlegrare the esuerione 
of motion numericaliy cver 46 long a tims interval as ig required, 


akon of the position af 


E 


and thus obtain o oumericoi repraceant 
the vehicle throughout the specified Cims interval. The limita- 


tions on this method Salli lato two classes. 
Tho first limitation is our imeccurete krowledge o£ the 
nysicei constants of tha system, Although it de beliewed that 
& 


e op de zt c- ^C s un - H % í ^ ` 
cur egtinebtes oZ the curthemcon mega regio may be in erres by 


as such ag 0.3 per cent, cut preseuk tings reeket experiments 
should significastiy esvones sur knowledge In ERIS ares. in 
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A a . s M H LN. "u > ej P An on Pe A T us + à A 
any case, tale Le nekt a basie Limitation of tha investigation 


of stebility. EZ wa can vhtain the answers te She stability 


thes 

questions baved on a model of tha earth-meen system where watt 
tà 

i 


la exzactiy 0.0123, we wiii have kaken a lona step toward 


el 
it 


obteining answers For the actual systen. later, with more 


precise knowledge of the mags otio, wo feel shat the reaulte 
do mot apply to the actual system, the procedures for compu- 


tation will stiil sppi: 


the major limitation on determining s precision orbit ia 
the errors inherent in the method of computation. The effect 
of the choice of a particular methet of computetion is dig- 
eussod by Baker end Gthersa (6). Basically, these errore may 
be classified ac truncation Gzror snd rounding error. Trunce- 
tion error rocaite from the use of s finire interval cover which 


p 


the valuc of the lategrand is epprseximated. It can be reduced 


tts 


$ 

by decreasing tho seize of the interval, end thus more ciosely 
approximating che vaiue of the integrand over the range of 
integration. Rounding orror is due to the limitation of compu- 
tation co a finite mber of disita, end the error ía cach step 
of the computetion io propagates throughout the entire problem. 
The reguiting coundínz evroz Shen tuexeases with the number cf 
SCHADE. 

These orrors can bs raduced by intelilgcaz selection of 
methods of computation. For example, For an orbit of this type, 
it is feit ae cusmpucation of pevturnations frau s To$erence 


2 £ rei £ n e e 
orbit weuld be mere afficsient than dixsask integration of the 
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ee 


Ge ce meee taS, sj mr oe e g: Dr ^ A 
acuatious. Howeyar, the cambinecion of truacation error and 


roundiug erzor Still esn not be eliminated. For aay precision 


{ 


orbit computetion, thera ie a practical limig to the duration 
o£ tho orbit which we can compute without a loss of validity 
ron the cumulativa correr. Thus, for a general lunar probe, 
zegardiess o£ rofinenents in the methed ef computation, we 
can specify where it will be on icta firs: pags near the moon, 
we cam entimate whexe it will be a fer months later, but we 


have no idea where it will be several years Later 


latuitively, it i8 feit that wa cen do better than thie 
with a& periodic orbit such ag ouga, Ik Seems that, 1f the 
orbit is periodic, we scan compute e precision orbit to complete 
one revolution about kke Trojan point, end the orbit fo 5tablo 
af it pesses closer te The Troien point the second Cine around, 
But we hove neglected the mation ef the vehicle in two modes 
with incomusnanuxvabia periods. & hypotheticel vahicie moving 
in the linearised syeten will or will nek meet the above eri- 


¿tea isdtial conditions, Bus wa know that 


pues 


terion according to 
the orbit of this vehicle is stable. Ie remains wichin upper 
and lover bounda determined by the initial conditions, and wil. 


pass arbitrarily close to any peint in thie region Lf we wait 


2008 enough. 
The next thoughts is that we should then compare tha computed 
precision orbit wita che colu&sios of the linearised equations. 


But again this compovisun given mo conelusive evidence of sta- 


tp 
n 


bility. As env comparison potat, the motion of the wehicl 


ES 








the actual 5ystew will hava altaren the proportions of rele- 
tive energy between the two modes. Wo may find thas tbe ac- 
tual rbit i9 closor to che Troiam point, but LE we compute a 
new linearized orbit using en imitisl conditions the actual 
displacement and velocity at the comparison point, we may find 
that the predicted upper bound of subseauent displacement from 
the Trojan point is greater than that predicted by the first 


linearized orbit. 


We can, of course calculate the predicted uppar bound and 
continue, wetehing chis upper bound fer en increase ox deezaass, 
But now wo exe back to tho problem o£ noc being able to wakeh 
it long enough, because ef the cumilarive ezzox o£ the numeriaal 
calculations. it is thie cumvlahiva error thet may mask the 
Stability effect that we wish to determine, eapecialiy LE the 
instebility of the orbit is ameli. We am working im a regien 
o£ cisliumar apace roughly boundad by two sllipsesa dofined by 
the solucion of che Lincarised equations. Since we ses no sig- 
mifieant pattern in the distribution cf the residual forces {or 
mon-lineay effeara} im thia regien, we should cover the entíx 
vegion rather unifomly In eux caieculetions. Otherwise, we 


cannot ba ture ther the west Few revclutiens of the vehicle 


es 


gbout Che Trojan point will not pasa through points whore the 


nog-iinear offoctso sre entirely diksereut. 


length of ¢ime nacesenry for an orbit to pass wichio 


& apecified distance e, 0, oe every point im its permitted reiiex 
varies with tho diszeibution of the sotion behween the two modes 





o£ pericdie motion. IS the motion o£ the vehicle 1a primarily 
confined to cre mods with the oiher ef empiisude no greeter 
than é then the necessary time Ls es longer them the period, 
of the primary made. TE tho amplitude of the secondary mode 


18 large, tha time way be prohibitive due te cumlative error. 


This leads us to consider investigation of the special 
cases where, besod on tha linearized model ef the system, Che 
anitiel conditions axe Such aa to excite onir one mode of 
vibration. The trouble with this is chat we should expect 
the non-linearities of the actual system to canae transfer 
of energy to the other mode» with a cumuletívo effect until we 


re back with the general case again. 
it seems possible thet we may obtain usoful information 
by constreining the orbit te the simple ellipse, and inte- 


grating the effact of the residual force over ome revelution 


about the Trojan point. Intultively, ib is felt ches a net 
Amwarü Korea would be an indication o£ stebility while a net. 


oumwarcd force would indicate Imetebilisy. If this were done 
for orbits of various sizes, we would then have an indication 
o£ stability versus average distance from the Trojan pointi, 
This investigation has disclosed mo mathematical justification 
for asserting tho meaningfuinesa of this “stability indication’. 


However, it is felt that the diecovery of a valid stability 


eriterion wiil requíre calculations of thie type. 
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CHAPTER 8 


LINEARIZATION OF THE EQUATIONS OF MOTION 
AT POINTS DISPLACED FROM THE TROJAN POINT 





6.1 Linearization of the Equations of Motion at the Positaon 
of the Vehicle 





In a digital computation of the orbit of a vehicle about 
a Trojan point, it seems most efficient to use a reference 
orbit and compute deviations from this orbit. For a reference 
orbit, the first choice that comes to mind is the solution of 
the linearized equations ae developed in Chapter 5. However, 
for a small segment of the orbit, we should get a better fit- 
ting reference orbit by linearizing tbe equations of motion 


at a point ou the orbit, rather than at the Trojan point. 


Linearizing equations 4.3 and 4.4 at the point (a,b) 


rather than at the point (0,0), the Taylor expansions are 


des 32. 
Su. a È (a,b) + sy (a,b) {x-a) + 2 A (rob) 
YX 


Ox Ox 

(8.1) 
dU a 22 (a,b) + au (a,b) (x-a) + (a,b) (y-b) 
dy = dxdy 5 l 


(8.2) 
4.8 





For convenience in writing the equations, let R4 end R, 
be the distances from the point (a,b) to the earth and moon 


respectively. Then from Figure 4.1 


RI = (b #3) + (a + e b? £p a* tatil (8.3) 


go b(a: ien abt k ^ -a41 (8.4) 


Equations 8.1 and 3.2 can be written a8 


a = k, +k, (x-a) + k (y-b) (8.5) 
du : a ü i» 

E ek, + kg (zea) + k, (y-b) (8.6) 
ae 


where, from equations 4.17 through 4.19, with substitution 


of equations 8.3 and 3.4 


k, = m (a,b) — -Cath -A (ai) R17 -6(a-) R5" (8.7) 
2 E - 3 2 -5 ec5 
ky = £y (a,b) - 1- ARA" - 68, + 3A (aty) BR,  3€(a-4) Ro 
(8.8) 
2 fa  -5 : -5 
ky = Lu (a,b) 7 3A(at5) (b+ 2) Ry + 36€ (a-2) 3) Ro 
KOY e 2 
(8.9) 
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pr Year Barr Da .« y (8.10) 
^ dy "de oe q? 


-2 -5 
ae 3 -3 E d, 

ks = 32 (a,b) e» l1- AB, - e, * 3A * 2) R} + 25 (b +3 Ro 
(8.21) 


ín linearized form, egostions 5,3 end 4.4 cam be written 


É -2 y» (1, > de, - blo) + kox + hoy (8.12) 
y 2x9 (k,-ek,- bk) + kys t koy (8.13) 


These equations can be written in differential operator form 
as 
(p^ - k) x - (p k,) y e kj - ak, - bi, 
(2p - kg) æ+ (p? - kj) y 9 i, - ak, - bis, 


The cheracteristic equation of this gysten is 


9 
(p° = kg) (p” = kg) + (2p + kg) (2p = Rg) =O 
GE 
F a 
pP t (4 - kg - g) P * (ks - ks en O 


and the roots of the characteristic oquatíou are given by 


E o 


e neis i ea = 
9 ENTE ARES ACRE me wenn 
& & 





ng 


= 


Jg Clo cae EPL, sD 


2 
E) a A 
ic- (Kalk 3 mka) 
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A stable solution is obtained if 
2 
(5-k, kk.) 
GD o 2 - 
in Ck + 5 7 0 
Substituting for the k's from equations 8.7 through 8.11, 
the stability criterion is 


2 wi, e) 9e A (p 413? 
yg a7 
ee 

(8.14) 


Equations 8.12 and 8.13 can be solved to obtain a solu- 
tion in the form of constants plus pure 8inusoids if the 
criterion of Eq. 8.14 is met. If the criterion is not met, 


the solution contains exponentially increasing sinusoids. 


Whether or not any significance can be attached to this 
criterion is not immediately apparent. If the solution resulting 
fron a linearization of the equations about a point displaced 
£rom the Trojan point is to be considered better than the golu- 
tion of Chapter 5, it can be used only over a small segment of 
the orbit, and then a linearization in the vicinity of a new 
point must be accomplished. We cannot then categorically 
state that the stability of a specific periedic solution is 
significant if we consider only a email fraction of one period 


o£ the solution. 


On the other hand, it is felt that, if the vehicle passes 


through a region where the solution of its linearized equaticns 
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is unstabie, there should be a tendency towerd stability o 
the complete orbit. If seems that each time the vehicle 
passes through this region, there would be & ewselative ten- 


dency for it to be pushed further away from the Trojan point. 


Time was not available for a survey of the region in the 
vicinity of the Trojem point to determine where the criterion 
of Eq. 8.14 is met.  Roughly it seems that it will uot be met 
for either x ox y graster then about C.1. In any esse, we 
Should £irat have a better justification cf the significance 
o£ such a survey, Although a time Limit on the present study 
preciuded further investigation along these lines, it ie felt 


that there may be something useful iu this area. 


un 
M 





CARRER 9 


FORGED TIS WW, RIERA obi 


OSN METHODS OF ORBITAL DES ERMIMAT Ti 
v$4 vato? R7) 4253 DOCET OS SAW FE 9 & AS EIA PDB IE LL IS ET DE ET 


At thie poing, the ignorance cf the suther of the Íntri- 
cacios of classical eolestial mechanica ia probably already 
apparent. The literature has been searched, but not thoroughly 
digested. If the meshed by which cur problem will be solved 
remains hidden in thie field, it cen moat probably be found in 
the work of Aurei Vintuez (7). 


: The recent Aussien worls on the subject: nag not bann Eho- 
roughly exeminod by the euthor, partiy because of the non-auailse- 
bility of teengletions, The work of Mermen (8) is beiafly dia- 


euaBgod by Leiaaeie (9). This morhod seca €e be alee described 


os 


by an abstract of a paper by Schigolev (10): 


"The motion o£ am asteroid is caleuleted acearding to 
approximation fovmulas..... The aspronimation values desizved 
are obtained by averaging. The gverege method cf the euther, 


however, Gifiers fran the woli-knowa mekhed cf celestish 


~e, $ Aw go 9. a? sz Par ; ^ Na 2 Ort T. a y . a «' 
gechsoicg inesuoch as The author dees nek average Che power 
fimetion, but contain stlifuilys choses parto of the dif 

5 à 5 e e SV ^ 4 e e "E 4 e ; sn e r s 7 
tial equations sud thus obtaius solvable differential equatiíoas”. 


fw, D n Ta caua : A onte Qe 4 fa ÉS prat i : > a 5 "f Lo 
his method wasn isvestigeted far ansugh to ascertein that 


peviedie solutions were indeed obtained, but not far enough to 


judge the walidity of these solutiona. 
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CHAPTER i0 


O EE CER 


IG Ad M a YS 


wb Y 1 


* 


^ 


4a stated ln Chapters 1i, ve do not have gpenific enewors 


to the problem under aoneideration. I£ we are aked "How 


a 


* ele re e Y, > e E i g Ads 
the orbital stability of ea peoposeó Trojen point satellita 


iu 


bo prediotod We csa fuegos Cho follewing arpreachy 


i. Evaluate the expected guidanee lnasaueacy Ín texme 


opf velocity sud dierlacement: error ab the Seoiam POLLE., Seller 


ehe I paiak es Surzested ia Seotion 5.0 to minimise the 


size of the ree auitiog orbit about tiie Prota mi point, based on 


‘4 


Ehe colusion of the linearised equations. 
£5 on Sema e my rd "em pg. “gi E ne en „= un 
4. Compute the orbit acound the Trojan roine, using as 


lneremental referente orbíta, golutíton 


i 
Er 
al 
ST 
ee 
i 
aT» 
Vw 


ecuscions of 
Chapter 8. The Lem ügth of Che inerements used in The compu- 
Cation will be picked to winimine the combination of truncatios 
error sud rounding exror. Also an atteupt should be made to 
vediet how tho sine of thls erxor wLil lucrenso wich time. 


* y A. sd Ea 2° to Qin 4 - os Loa AE as = a Eur e - 
LE Ehe comoutod oxbit divevrgosc Eran the * zojan noise before chis 


t El * N A o a e Aze Sek m. $ p ` 
computation exror beccues Sigodficane, tbe quention ism onswered 


es fo — «e^ y - Aa EN e ¡E nA 5 ^ 
But a correct xecsbhva suswer is much enniex Ga ebt@gan Ghar 
wor 2 en de c 7 T NA TE Kein ie cog 63 C^ 
& corzreot boutcivo auawer. IR the eoupated ochic df verges Laws 


the Trojan point omiy siightiv 


a 
e^ ^ CONV en, mE - QU 
wer CEN a 





g 


BE» a on: a? e 1 £. . E au, RR Let ^s x$ > ES 
thia is dus to aa actual imetasiiity ok Cho OXOLC, GE a result 


ck the acoumulatad commutation error. 


3. Further work on tha orbital stability of Trojan 
/ ) $8£ellites should probably be directed teward the problem of 


identifying and zedueing the effect of computational errors. 


Aliso it is felt that further investigation should be made 
o£ the significance of the region In which stable solutions 
of the linearized equations axe obtained, as diseussed in 


Chapter &, 





APPENDIX À 
ANTENI IIIA AA ERA CC e REA 


NUMERICAL SOLUTION OF SER LIUESEIZED DOUATTOIB OY MOTION 


The foiloving calcnlations were mado to pexmit cgicu- 
tation of a golutica of the iinesrlized eguarions diracti: 
from specified initizd conditions. For nmumeaivt.ecal werk, s 
value of © of euraathy 0.031213 mas used. Aithrugh, the 
physical quantity $a Sha eartlı-moon system Le only hvown 
to within about 0.3 por cent, the caloula ion), were cevriod 
to ton significent Figures Lo permite us? Of the resulte la 
Subsequent numerical Integrniien compukiitions, 1£ desired. 
Tbe results are then conmsintent with . linearized model 


ef the system with the macd ratio epelifiaed by € w» 0.01213 


fo @ i. i 2 0,98787 


11-276 (1-6) e 0Q,82:65735232 





2 c » y z 
po 7 & (190.8226735232) w 0,0807632386 


Pa = 0. 9545872205 
Po = 0.2979318002 


ib 


(1-28) = 1,207523462 


Substituting initial eonditions = (0), y (0). % (0) 


and y (0) into equacions 3.1 and 5.2 


& (0) = pzAy + 190, (A. 3) 
y (0) = pA. + poto (4.4) 


Gombiaing equations 5.3 end 5.6 we cau geb A, iu 


Tarma of 8. aud Ba 


Un 
ved 





t 


dp, v An. (1-26) B, s (pj + 7) Be (4.5) 
Similarly fer equetions 5.5 and 5.0 
29940, > 3Y3 (1-28) D, + (pf + 2) Dy (4.6) 
& és 
From eg ations Aael md 4.2 
D, cx (0! - By (A.7) 
D, = y 0) = By (4.8) 
Jubecitusing oquetiono &.7 end A.B into Eg. A, 


€ T | i | 
PPa a E (1-25) |x (0) - 5 X (va y 2 zos] 
(4.9) 


Multiplyio; Eq. |.3 by 2 aad substituting from equations A.S 
and 4.9 


è Z5. 
ae 0) 223 Cane) By + 5(0)-B, [+ (nye yaybtnd + ©) [¢0)-8, | 


(A, 10) 





2x (0) © aia (1-28) 2 (0) +3 = Bg) Byrtoge 2960) 
» a. (4.11) 
X3. B & 
By E aq9)- opem (0) *-$—- g & (9) 
PT Po By 7 Pa Up UR 
(A.12) 


Subatitkuring mmerleal values 


By = -1,5811113%4 wid) -2.843972877 y(01 4 ?.431689219 x (0) 


* 


(A13) 


Gonbining equations 5.3 and 5-4 te eliminate A, instead 


of Ag, We can get i, im terme of E, end By 


= TIE ire | 
ep 155 - > a + A By - UA € Loge } Ba ¿A # 16) 
Similariy fsen equakions 3.3 and 5.6 


ada | 
Pot a e (pi 4 3j 24 P dat: (1023) By ¿£0135) 
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Substituting equations A.7 and A.Y iuto Ix]. A.15 
> G Mu 
ge mr e al a ms sef = vx 5 us po TA PE i AGIT 
Es ( Yo nU ES p, | , (1-2€) | y(u) Ba | 
(A.16) 


Multiplying Eq. A. by 2 and substituting from equations 
A.16 and A.16 


2y (0) = - (pj + $ By > 38 -2€) DE y) - 5 | 


y (8.17) 
AE | — a 
2y (u) = - 38 0-20) yO) - ep By ~ (3 + 2) sO) 
(A.18) 
»> +3 ifi (1-2€) 2 
B æ San - : = Y ee dtp RE 
B4 Aare x(u) > 5 yn 5 y(u) 
u 72 Pi 7 P2 Zu 
(A.19) 


Substituting numerical values 


B, = -1.019805762 x(0) -1.56111154645 y(U) - 2.431689219 y(u) 


(A.20) 


Gu 





“se & x <+ . Md 
fy can now be found from Ed. A. 


7 
T 
A3 2 
i 


A 


» 2 
f 22" id 3 y 
A. m uu Alee) — B4 2. bea T By (4.21) 
2p4 epa : 
Bg = 0. 6659118009 By + 1,055213473 Bo (A.22) 


Substituting from equations A.20 and À.13 


AQ 3.22885432380 x(0) - 3.731583091 y(O) 


* 


44.026523771 x(0) - 1.615427169 y(O) 
(A.23) 


Similerly from Eq. A.14 
"a 3 343 ~ 
a- ow f > i- E 
Ag E By | P. Bo 


(A.24) 


Ay = -0.8701335645 B4 - 0.6639118009 Bo 
(4.28) 


Substituting from equations A.20 end A.13 





A, =» 1,910529363 x(0) + 3.228056338 y£0) 


-1.614427169 x (0) ^ 2.115894408 y (0) (4.26) 


Substituting Eq. A.20 into Eq. 4.7 
D, » 2.019805762 x(0) + 1.543111544 y(0) + 2.431689219 y (0) 


(4,27) 
Substituting Eq. A.13 into Eq. 4.5% 


D, 9 1.541111544 x (0) b 3.843572077 y (0) - 2.521699219 x (0) 











(4,23) 
C, can be found from Eg. 4,6 
33 Ei 2 
c» 4 D. r EN e Dp 
Cy A bo 5 2 
p (4.29) 
C, * 2.127205440 D, 4 3.925000296 n, £A.30) 


* 
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I 


substituting from equations A.2/ and A.23 


C, 10.345£0507 x(0) # 13.366236/6 ye) 


-9 544390904 #(0) + 5.1727025: 


Similarly from Eq. A.1i5 


2, 3 ined 

e B cum {1-26} 
20 apy Pa a Pe 
C, e - 1.5076453965 D, - 2.127205440 D, 


Substituting from cquations A.27 and 4.29 


Cy = = 6.121428251 x(u) - 10.34540597 y(u) 


+ 5.172702535 x (u) - 3.422952654 y(u) 


35 y(u) 


(A. 31) 


(4.32) 


(A.33) 


(A.3%) 





2 


Given specified LoítiaelL canditioss, the constaate 


a 


of integration can now be caleuleked fron equations 4,13, 
A.Z0, ÀA.23, À.25, A.27, 28. Avril and A.M. These equa- 
tieas were crosa checkad nunarically by ether formulas and 


found to be conaistent. 


Tho auberítution ef muabers into this equetion talla 
ug little about whel sort of an achte we would get with 
slightly cizlexent Initial conditions, as does the nathod 
of Secticun 5.5.  Houcwer, Jor saleulatier af a speeifie 


L ^ 


cONDL6, Cho mzerioah work e mush simpler- ££ this solus 


tion wore required al a Kutetense orbit for a digital 


ecmputationm o£ the orbit, the aboye sqoatiovro would be used. 





UME. LOCAL COMP) TAT ALLON í FOTE RESIDUAL JOCOBLAN WUNGILON 


To determing soma numerical values of u and u', 
erations 6.5 end 4.20 wexe transformed to polar coordi- 


mates by Kin equacicnas 


The eauatioas then bacaue 


y”. € Deos 94 y k (p^os ga“ Oa\3 ania det co 0 + peos Ba) 


ed 


3 ar m a 
GA punto, 30 sinds P e cog O4 E cos É 4- i) 


2. 
e AN 3 = 
tí p oin As GG Q gin’ + p? ces” ~p aaa Or i) 
3 
(Beh) 
qoare» 
(2 UE ' ` 2 9 2 2 
13 zu 3 el 3 y & ca e Y Buc { mee ya Sin b con Ú + és AS e 
2 8 a 
(2.2) 
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Using che idontitios 


sin“ Q cog em i 


Y Ban Ó + 4 eos Da sin 60° sin Ús cos 60% cos 9. 


cos ( G 60°) 


gin E 5 cog Ü s. E25 (E 60") 


Ro doo 


The equations Kor u and wm! become 


3 


ur -© 9.0504 à k Jr * 2 ("cog (0-609) & 1 


awe 


ul * 29co8 (U-60%) + 1] +e y ^ «2 8 eost 07e 
eun "e 
(8.3) 


; P À 3 * m 4 » d 2 
um 34 2 Oo” + WM (1-26) 0 ?ain & G.. 2 9? con 28 
3 dy e G eh 3 


vie 


(B.5) 


Using e value for G of O.Ol213, Ehe quenzizies u and u' 


wore enlcuietco0 for 0.2 imeramesta of feroz O te 0.5 mad fer 





30% incremento of 6 , The resulte axo listed in Table Bal. 
s + ? 
Tho zegSidusi function u = u- c! at these polate £e ohowa 


in Figs B, ay 


p 
Vul 





teo? 
457 


8 224 e efto 
ise” »-82 p? 
08 e 204 - 
o 9/ ob, en {SS 
of 
237 n ! tS 8 
o 4 ed (73.29 safe 
© i i ia 
$80 ee6G29 1-322 0-35 «033 «$ 0.0 o 4 a2? o42 o&zo RA D? is 
3 E 
0 a cmi 
025 ED eS y o 6 
9 go o 
997 sang 29 
? (e$ s d 
E ats” ^ 279 2713] o ES = 
210 $0 330 
3 ¿(OSO 310 
o AES 
o 295) a -603 
et = c 9t 406° 
215595 


£0? 
oe G/F 


Distanees from the Trojan voir im 0.1 inczesencs 


v 


T * 
indicated values ere . x ant 


d cO 


Fig. RA, 
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Approximate Values of the Residual Jaccbien Function 





TABLE B.1 


XIABULATION OE JACOBIAN FUNCTIONS 


ETS "i i 
| ju | E 












2.5041 1.5033 
1.5177 | 1.5150 
1,5429 | 1.5339 

| 1.5825 | 1.5600 
5 | 1.6297 | 1.5938 





Co ee oc? trip IRA a pL asl 


1209| 1.5043 | 1.5039 % 2 
1.2 | 1.5183 | 1.5156 § 
| 1.5441 | 1.5350 
1.4 | 1.5946 | 1.5622 E 
| 1.6329 | 1.5972 | 
150%.1 | 1.5002 | 1.5001 


1.5010 | 1.5006 
1,5059 i. 5012 
1.5103 ! 1.5022 
1.5222 | 1.5034 


Iris eng ws mes 
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A MA | | 2 


HA nado Mg m po Pp 
i 


6 





= EJ 


n 
180° | 
to 


11 
.2 


+. 


E 
| wd 
j of 
| 2 


au: m 
1.503311.5038 


1.5112:1.5150 
1.5203 1.3335 
1,5273 1.5600 


| 1.5310 1.5939 


| 1.5114/1.5111 
| 2.3470] 1.9465 
1.60801 1.6000 
1.689441 1.8778 
| 1.802311.7778 
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